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I . Introduction 


Among  the  various  alternatives  which  have  been  proposed  for  the  propor- 
tional navigation  guidance  law  in  missile  intercept  problems,  Deys^  [1]  applied 
optimal  stochastic  control  theory.  Ho  [2]  and  Rajan  [3]  applied  the  theory  of 
differential  games,  and  Slater  C43  and  Nazaroff  C53  took  the  perturbation  point 
of  view  in  applying  linear  optimal  control  theory.  In  this  paper  a determinis- 
tic approach  to  a two  dimensional  missile  intercept  problem  is  considered  under 
the  assximptions  that  the  pxirsuer  possesses  thrust  modulation  capabilities,  in 
addition  to  the  usual  thrust  vectoring,  and  that  continuous  measurements  of  the 
relative  range  coordinates  are  available.  The  thrust  modulation  capability 
facilitates  the  formulation  of  the  control  problem  as  a commutative  bilinear 
system  allowing  for  the  application  of  some  results  recently  obtained  in  [63. 
The  philosophy  for  deriving  the  overall  control  law  is  basically  the  same  as 
that  used  by  the  authors  in  [73  for  the  terminal  control  of  a gliding  parachute 
system  in  a nonuniform  wind,  viz.  independent  and  easily  computable  solutions 


to  the  control  and  estimation  problems,  relative  to  a time  interval  t^  < t « 


are  combined  to  define  a step-by-step  control -estimation  sequence,  i = 0,1,2»», 
which  constitutes  the  closed  loop  control  law. 

Following  a statement  of  the  problem  in  Section  II , a least  squares  esti- 
mation scheme  is  developed  in  Section  III  to  estimate  the  target  speed  and  rela- 
tive heading.  A closed-form  solution  to  a minimum  control  energy  problem  with 
terminal  constraint  is  obtained  in  Sectibn  IV,  assuming  the  parameters  for  the 
target  speed  and  heading  are  known.  A modification  is  considered  in  Section  V 
to  include  actuator  dynamics  for  the  thrust  vectoring^ehg4ne . Simulation  re- 
sults are  presented  in  Section  VI  which  combine  the  estimation''p^ameters  from 
Section  -III  with  the  controls  obtained  in  IV  and  V for  the  closed  lM^-(feedback ) 
control  law. 
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II.  Problem  Statement 

It  Is  assumed  for  a hl^  speed  pursuing  missile  and  short  Initial  range 
that  the  maneuvering  of  the  vehicles  can  be  restricted  to  a two  dimensional 
plane  as  shown  in  Fig.  1 with  the  coordinates  fixed  in  the  missile.  Denoting 
the  angular  rates  of  the  missile  and  target  with  respect  to  a nonrotating  ref- 
erence frame  by  u^^  and  u^  respectively,  the  kinematic  equations  of  motions  are 
decribed  by 

" -’t  "s  ♦ *2”h 

ij  • cos  «3  - (1) 

*3  " “t  - "m 

speeds  of  the  target  and  missile  relative  to  air, 
Xj,  and  Xj  are  the  horizontal  and  vertical  distances  of  the  target  relative  to 
the  missile , and  Xg  is  the  xatlve  angle  between  the  headings  of  the  missile 
and  target  [5]. 

With  a specified  sequence  of  time  instants  t^,  i=0,l,2*»,  the  estimation 
problem  relative  to  the  i—  time  interval,  t^  < t < t^^^,  is  to  obtain  estimates 
of  the  parameters  (v^,u^,Xg(tj))  based  on  continuous  measurements  of  the  posi- 
—tion  coordinates  (Xj^(t),X2(t))  and  a knowledge  of  the  missile  line  speed  and 
turning  rate  (Vj,,Uj,)  on  The  control  problem  relative  to  the  i— 

time  instant,  t^^,  is  to  obtain  a reasonable  control  strategy  for  and 

t > tj,  such  that  an  intercept  occurs  at  some  fut\n?e  time  T > t^,  i.e. 
Xj^(T)  = Xjd)  = 0,  assuming  that  the  estimates  of  the  parameters  (v^,Uj,Xg(t£ )) 
are  exact  from  a previous  subinterval.  This  pair  of  estimation  and  control 
problems  is  solved  anew  on  subsequent  time  intervals,  t^^^^  < t < etc.,  in 

defining  the  closed  loop  control  law.  Precise  statements  and  solutions  to 
these  problems  are  given  in  the  following  sections. 
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III.  Least  Sqiiares  Estimation  of  (v^,u^,Xg) 

Let  denote  a typical  time  interval  over  which  continuous  measure- 

ments of  the  range  coordinates  (Xj^(t),X2(t))  are  assumed  given  together  with 
a knowledge  of  the  missile  turn  rate  and  line  speed  A deterministic 

least  squares  estimate  of  the  parameters  (v^,u^,Xg(t^))  is  considered  hei«  in 
order  to  obtain  a simple  closed  form  solution  to  this  aspect  of  the  problem, 
thereby  facilitating  an  easily  implemented  control  law  for  the  overall  missile 
intercept  problem  as  discussed  earlier. 

Denote  the  target  speed  v^  and  initial  heading  Xg(t^)  by  the  parameters  v 
and  "a",  i.e.  (Vj,Xg(t^))  = (v,a),  and  assiane  for  the  time  being  that  the  target 
turn  rate  is  a known  constant.^  A least  squares  estimate  of  (a,v)  results' 


upon  minimizing  the  functional: 


J(a,v)  = 


f ^ + V si 

I ^ - V CO 

•'t 


sin[a+U(t)]  - u^^-Ct)}-  dt 


V 

X2(t)  - V cos[a+U(t)]  + 


where  U(t)  is  defined  in  terms  of  the  target  and  missile  turn  rates,  and  Uj^, 


'1^  CV"h< 


s)]ds. 


A necessary  condition  for  the  minimization  of  J is  that  the  partial  derivatives 
of  J with  respect  to  "a"  and  v vanish: 


= 0 
3a 

(a*,v*) 


(a*,v*) 


^The  estimation  of  u^  will  be  considered  separately  at  the  end  of  this  section. 
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Observing  that  J is  quadratic  in  v,  the  best  estimate  v*  can  be  uniquely  deter- 
mined in  terms  of  a*: 

^ ~ ■ Xj^(t^)8in  a*  + X2(t^)co8  a* 

1 o ' » 

- X2(tj^)co8Ca*+U(tj^)]| 

- [ ^ {uj.|x3^(t)cos[a*+U(t)]  + X2(t)sin[a*+U(t)]  + cos[a*+U(t)]|  dt. 


Similarly, 


whex^ 


0 = A cos  a*  + B sin  a* 


A = cos  U(tj^)  + XgCtj^)  sin  U(tj^)  - Xj^(t^) 

+ jcvj,  + lycj^(t)]  sin  U(t)  - ®o® 


B = sin  U(tj^)  + ®°®  U(tj^)  - XjCt^) 


I ^ jcvjij  + u,j,Xj^(t)]cos  U(t)  + sin  U(t)|  dt. 


ri?om  Equation  (6),  we  obtain 

m = 0,±1,...  (7) 

This  is  equivalent  to; 

B cos  a*  > 0.  (8) 

Following  the  solution  for  a*  from  (7)  and  (8),  v*  is  obtained  from  (5). 


a*  = 2iinr  + tan”  ^ 


2 2 

in  which  m is  chosen  so  that  9 J/3a 


(a* ,v*) 


The  above  solution  for  estinatlng  assumes  a knowledge  of  Uj.. 

Hinimizlng  (2)  over  u^  will  not  lead  to  a simple  closed  form  solution  and, 
therefore,  another  approach  is  necessary  to  estimate  u^.  A new  approach  to 
parameter  estimation  problems  has  recently  been  developed  in  [8]  which  can  be 
directly  applied  to  this  problem.  In  applying  this  technique  it  is  first  neces- 
sary to  rewrite  the  underlying  differential  equations  (1)  in  terms  of  a dif- 
ferential operator  equation  of  the  generic  form 


0 = P(D)w(t)  - Q(D)V(t)f(e)  , t < t < t,  (9) 

O X 

where  P(D)  and  Q(D)  are  polynomial  matrices  in  the  differential  operator  D * ^ 

dt 

given  by 


" iL  ‘“-i  “i 


Q{D)=  I 

i-O 


and  (w(t),V(t))  are  vector  and  matrix  valued  functions,  respectively,  of  the 
observed  data  on  [t^,tjj^].  The  vector  valued  function  f(e)  in  (9)  is  a given 
single  valued  continuously  differentiable  function  of  the  parameter  vector  9 
which  is  to  be  estimated  based  on  the  given  data  over  Ct^,t^].  For  the  missile 
Intercept  model  (1),  the  requisite  equation  can  be  obtained  by  differentiating 
the  first  equation  in  (1),  substituting  from  the  second  and  third  equations  in 
(1)  to  eliminate  v^  cos  x^  and  x^,  and  rearranging  the  resulting  expression  to 
obtain  the  following  differential  operator  equation: 


0 « CD^,D,1] 


-Xj^(t) 

X2(t)[Uj^(t)  - 1] 

- CD,1] 

Ujj(t)Cx^(t)xij^(t)  + Vj^(t)3 

Xl(t)Uj,(t)  + Vj^(t) 

• « 

(10) 


This  model  is  in  the  form  of  (9)  with  P(D)  = Row[D  Q(D)  * RowCD,13  and 


rvBj  = B » \i^  - tne  parameter  to  ne  estimated  given  the  data  on 
Ct^»tj^].  Applying  the  theory  developed  in  [8],  it  is  then  possible  to  obtain 
an  explicit  expression  for  the  estimate  without  estimating  the  unknown  ini 
tial  conditions  (Xj(t  ),x„(t  )).  Since  this  theory  was  developed  after  the 


current  work  of  this  paper  was  completed,  these  expressions  will  not  be  written 
down  here  and  the  simulations  in  Section  VI  will  only  incorporate  the  estimates 
for  (v_,x-). 


IV . Minimum  Energy  Control  of  the  Intercept  System 


Let  t^  denote  an  arbitrary  initial  time  with  estimates  of  the  parameters 
(v^,u^,Xg(t^))  assumed  given  from  a previous  time  interval.  Rather  than 
consider  Vj^(t)  and  u^^Ct)  as  independent  control  variables,  a proportionality 
relation 


is  postulated  with  the  proportionality  parameter  Y to  be  determined  from  the 
boundary  conditions.  The  reason  for  this  postxilate  is  based  on  expediency  in  ob 
taining  a closed  form  solution  to  the  following  optimal  control  problem: 


Minimize  J(u) 


subject  to  the  intercept  condition 


for  some  finite  terminal  time  T > t^  (a  free  time  formulation),  where  the  nota' 
tion  u s U|^  is  used  for  simplicity. 

The  system  (1)  can  be  transformed  into  a commutative  bilinear  system  by 
introducing  the  following  aiixiliary  states: 


That  is,  with  these  additi^hT ^slgtes  (1)  can  be  written 


It  can  be  readily  verified  that  AB  = BA  so  that  (14)  is  a "commutative"  bilinear 
system.  The  implication  of  this  fact  is  that  the  optimal  control  solution  to 
(12)  and  (13),  if  it  exists,  is  simply  a constant  function  as  has  been  shown  in 
[6],  Theorem  3.  The  existence  is  contingent  upon  the  terminal  condition  (13) 
being  a reachable  point  for  the  given  set  of  initial  conditions.  In  order  to  ex- 
amine  this  possibility,  assume  that  v^  and  u^  are  constants  so  that  (14)  can  be 
Integrated  explicitly  as  follows: 


u(8)ds 


, x(t  ) = 

0 

sinx.(t  ) 

- 1 

3 o 

COSX-(t  ) 

0 0 

1 . 

Xj^(t) 

Xj^(t^)cos  j u(s)ds  + X2(t^)8in 
^o 

j U(8)d8 

^t 

o 

ft 

ft 

-x,(t  )sin  I'  u(s)ds  + x-(t  )cos 
' 1 o J.  2 o 

j u(s)ds 

O ‘o 


(17) 


-sin[x,(t  ) - 

d O 

•t 

uds  + u^(C-t^)] 

0 

ft 

<i5+ Y 

cos 

■t 

uds  -1 
t 

o 

ft 

cos[x,(t  ) - 

3 o ^ J 

uds  + u_(C-t  )] 
t * o J 

o 

-sin 

J 

uds 

t • 

o 

We  shall  resolve  the  terminal  constraint  problem  by  considering  the  inter- 
cept angle  as  a peu?ameter,  then  incorporate  this  solution  with  the  minimxnn  energy 
problem.  Consideration  should  be  given  to  two  separate  cases  in  which  is 
zero  and  non-zero,  respectively. 


Non -zero  Angular  Maneuver  of  the  Target 
The  terminal  constraint  (13)  on 


(Uj.^  0) 

and  X-  requires  (for  some  T > t ); 
^ o 


CXl(tp)+Y]cos[u^(T-to)+X3(t^)-6]  + X2(t^)sin[u^(T-t^)+Xg(t^)-e]-Y 
,-[Xi(to)+Y]sin[uT(T-to)+X3(to)-6]  + X2(t^)cos[u^(T-t^)+X3(t^)-3] 


cos  6 - cos[u^(T-t^)-e] 
sin[u^(T-t^)-e]  + sin 

T 


(18) 


where  6 * ~ \ ^ defined  as  the  intercept 


angle. 


Therefore,  the  terminal  constraint  problem  has  been  reduced  to  solving  a 


pair  of  transcendental  equations,  (18),  for  an  appropriate  set  (y,8,T).  A solu- 
tion often  exists  for  this  case  in  which  the  number  of  unknowns  exceeds  the  nxun- 
ber  of  equations.  Fixjm  Equation  (18)  we  obtain 


V t^'V-»->»vk^  *a>-i:t,W'V-a.'^,p^-*^-»'xj>.'-yvW,»^'rtJ-— w4.^. 


cose  . cosXgCt^)  . ^x^(t^)  t ^x^(t^)co8X3(t^)+X2(t^)8inx3(t^)  - 

(19) 

Given  {x, (t  ),x.(t  ),x,(t  ),u_,v_},  appropriate  though  nonunique  values  can  easily 

XO  aO  OO  11 

be  determined  for  y and  6 such  that  (19)  is  satisfied.  After  y and  8 are  so  de- 
termined, the  intercept  time  T can  be  computed  from  (18)  as  follows: 


T = t + ^2kv  + tan“^ 
o u^  G 


(20) 


where 


F = Y[Y+Xj^(t^)]sin[B-X3(t^)]  + y cos[e-X3(t^)] 

v_  \ 

- ^ (ysin  B-[Y+x.(t  )]sin  X3(t^)  + X2(t^)  cos  X3(t^)} 

G = Y[Y+X3^(tQ)3cosCB-X3(t^)!l  -y  *2^^©^  sin[B-X3(t^)) 

V 

(ycos  B+[Y+Xj^(t^)]cos  X3('t^j)  + ~ ^ 


k = 0,±1,...  such  that  T > t^. 


In  siunmary,  if  a constant  but  non-zero  angular  maneuver  of  the  target  is 

assumed,  then  there  exists  a triple  (YtB,T)  satisfying  (19)  - (20)  which  solves 

the  initial  and  terminal  constraint  problem  (13)  - (15)  for  every  (x, (t  ),x.,(t  ), 

i.  o Ji  o 


x_(t  ))  e R . The 
o o 


'boxires^nding  proper  control  action  satisfies 
fT 

J u(8)ds  « *3(^0)  - B 


(21) 


Zero  Angular  Maneuver  of  the  Target  (Uj,=0) 


)]} 


In  a similar  manner,  the  tezmiinal  constraint  becomes 
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’[Yi-Xj^(t^)Jcos[B-x^(t^)J-x^(t^)sin[B-x^(to)J-v^(T-to)8inB  - y 

0 r 

[Y+Xj^(t^)Jsin[e-x^(t^)J+x^(t^)cos[6-x^(t^)J+v^(T-t^)cos  $ 
where  B is  as  defined  in  (18).  This  equation  can  be  reduced  to 


cos  B = cos  XgCt^)  + i [x^(t^)cos  XgCt^)  + X2(t^)sin  XgCt^)].  (22) 
Finally,  the  intercept  time  T is  determined  by 

‘ '^o*  ^fY+x^(t^)3sin  X3(t^)  - X2(t^)cos  X3(t^)  - y sin  B).  (23) 

A feasible  T requires  T > t^,  i.e.,  the  term  in  the  bracket  must  be  positive. 

It  can  be  shown  (see  Section  3.2  in  [9])  that  this  is  ^e  only  for  those  ini- 
tial conditions  outside  the  region  E defined  by 

E - {(0,y,2)  e R ; either  y > 0,  z = (2k+l)ir,  or  y < 0,  z = 2kir;  k = 0,±1,...} 


Therefore,  if  a zero  angular  maneuver  of  the  target  is  assumed,  then  there 


exists  a triple  (•y,B,T)  which  solves  the  initied.  and  terminal  constraint  prob- 
lem (13)  - (15)  for  every  (»i<t^).X2(tQ)»X3(t^))  e \e.^ 

It  should  be  noticed  that  in  the  previous  analyses  the  control  function 
u(t)  has  been  eliminated  for  simplicity  of  computation.  However,  an  admissible 
control  which  steers  the  missile  to  the  target  at  T is  associated  with  the 

triple  (y,B,T)  through  Equation  (21).  Thus  the  set  U of  admissible  controls 

c 

is  specified  by: 


2 f’’ 

= {ucL  ([t^,T],R):  j u(s)ds  = X3(t^)-B+u^(T-t^)} 
where  B and  T are  determined  by  Equations  (19)  - (20),  or  (22)  - (23). 


(24) 


B denotes  complement  of  B in  A. 


I i 
1. 


f 
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With  the  set  of  admissible  controls  furnished  as  in  (24),  we  are 
ready  to  state  the  solution  to  the  minimum  energy  problem  (12).  The  following 
proposition  is  a direct  consequence  of  Theorem  3 in  [6]. 

Proposition  1; 

Given  the  system  (14)  - (15),  there  exists  an  optimal  control  u*  e U 

c 

which  minimizes  the  cost  (12)  subject  to  the  constraint  (13)  for  each  appropriate 
set  of  initial  conditions  ccntrol  is  given  by 


* X (t  )-6 

“ = “t  ^ 

O 

where  T and  B are  given  in  Equations  (19)  - (20),  or  (22)  - (23). 


(25) 


V.  Singularly  Pertiirbed  Problem 

A more  complex  system  model  is  considered  in  this  section  in  which  the 
missile  turn  rate  is  taken  as  the  output  of  a first-order  lag.  This  takes  into 
account  the  practical  situation  in  which  the  missile  turn  rate  is  furnished  by 
a DC-motor  having  first-order  actiator  d3mamics.  That  is. 


e u(t)  = -u(t)  + u 


£ > 0 


(26) 


where  u^  is  the  real  input.  In  the  limiting  case  where  e approaches  zero,  this 
consideration  is  generally  known  as  a ’singular  perturbation'  problem  [10]. 

The  cost  functional  in  this  case  becomes 


J(u  ) 
o 


-if  ! 


(s)ds 


(27) 


subject  to  the  constraint 


x, (T)  = x.(T)  s 0 for  some  T > t . 
i.  ^ o 


(28) 


.V 

) 


ri 


I 


By  defining  z = u,  the  system  equations  can  be  expressed  by; 


X a Ax  + Bxz 


ez  a -z  + u 


e > 0 


where  A and  B are  the  same  as  defined  in  (15). 

Because  it  was  shown  in  the  last  section  that  the  terminal  constraint 
problem  has  a solution  when  the  control  action  satisfies  (21),  or  in  terms  of  the 


new  state  z 


I z(s)ds  a XgCt^)  - B + u^(T-t^), 


the  terminal  constraint  problem  (28)  also  has  a solution  provided: 


l^T  rt  - — 


o o 


-]■ 


Thus  the  set  of  admissible  controls  for  the  problem  (27)  - (29)  is  the  col- 
lection of  inputs  u^  satisfying  (30). 

The  existence  of  an  optimal  control  u^  e which  minimizes  the  cost  (27) 

can  be  easily  established  and  by  the  Maximum  Principle  this  optimal  solution 
satisfies 


u^(t)  = q(t) 


p(t)  a -3H/8x  a -(A’+B'z)p 


cq(t)  a -aH/az  - -x'B*p+q 


« eq(T)  a 0. 


(32) 


In  the  latter  expression,  (p,q)  are  costates  corresponding  to  (x,*)  and  prime 
^dei^es^h^iBayx  transpose  operation.  It  can  be  easily  checked  that 
^^x'B'p)  s 0;  hence  x'(t)B'p(t)  = k,  a constant  to  be  detexihlHed  fay  the  beund»? 
ary  cmditions. 

Substituting  k into  (32),  q can  be  solved: 


q(t,c) 


(33) 


Then  from  Equations  (30)  - (31),  k is  given  by: 


T-t 


X3(t^)  -B  + u^(T-t^)  + ez(t^)(e  ® -1) 


(34) 


where 


L - c 


+ T-t 


These  results  are  summarized  in  the  next  proposition. 

Proposition  2: 

Given  the  system  (31),  there  exists  an  optimal  control  u e U which  mini- 

o c 

mizes  the  cost  (27)  subject  to  the  constraint  (28)  for  each  appropriate  set  of 
initial  conditions.  This  control  is  given  by  (33)  and  (34)  where  g and  T are 
specified  in  Proposition  1. 

It  should  be  noticed  that  the  singular  perturbation  comes  into  the  prob- 
lem (29)  as  e appraches  zero.  This  is  clearly  seen  from  the  expression  of  the 
optimal  control  (33)  - (34),  i.e. 


-14- 


11m  k(e) 
e-K) 


o 


lim  u*(t,€)  = - 11m  k(c) 
e-K)  cK) 


x,(t  )-6 

j o 


< t < T. 


which  Is  exactly  the  same  as  that  In  (25).  Therefore,  for  e sufficiently  small 
the  solution  to  the  optimal  control  problem  (27)  associated  with  a fourth-order 
system  (29)  can  be  approximated  arbitrarily  closely  by  the  solution  to  the  prob- 
lem (12)  associated  with  the  third-order  system  (1). 

In  fact,  not  only  the  reduction  of  system  order  Is  shown  here,  but  also 
an  explicit  solution  to  the  optimization  problem  of  the  quadratic  system^  Is 
derived.  This  provides  a great  deal  of  potential  to  Implement  such  a control 
law  In  practice  because  the  firs;,  -nxler  actuator  d]rnamlcs  have  been  Included. 
Actually,  this  result  can  be  generalized  to  Include  any  higher  order  actuator 

I 

dynamics  as  long  as  the  constancy  of  the  control  area  Is  sustained. 

VI.  Slmulatlotf^Results 

The  least  squares  estimation  scheme  of  Section  III  Is  combined  with  the 
optimal  control  law  of  Sections  IV  and  V to  form  a step-by-step  feedback  control 
of  the  missile  Intercept  system  for  selected  Initial  conditions  and  target  turn 
rates.  The  Initial  taz^et  turn  rates  are  assumed  to  be  0.0  rad/sec  and  0.15 
rad/sec  while  the  target  line  speed  is  lOOOft/sec.  In  addition,  two  sets  of 
sinusoidal  fluctuations  were  Imposed  upon  the  target  turn  rate  and  line  speed 
as  listed  in  Table  1 to  test  the  sensitivity  of  the  system. 

Note  that  (29)  Is  no  longer  a bilinear  system  as  defined  In  Section  IV;  Instead, 
It  is  sometimes  referred  to  as  a quadratic  system. 


A time  constant  e = 0.5  sec  was  assimed  for  the  first-  order  actuator  dynam' 
ics.  During  the  first  estimation  interval  of  duration  0.5  sec.,  the  control  ef- 


tervals  were  determined  by  Equations  (20)  or  (23)  (depending  on  whether 


of  control  led  to  an  interception.  Additional  estimation  and  control  intervals 


were  needed  for  the  perturbed  target  turn  rates  and  line  speed  cases  shown  in 
Figs.  4 and  6.  It  is  clear  that  in  the  ease  where  the  target  turo  rate  and 
line  speed  are  constant  (Figs.  2,3  and  5),  the  estimation  of  the  target  speed 
and  the  initial  relative  heading  is  exact  and  an  ideal  intercept  is  achieved 
in  one  interval  as  expected.  For  the  case  in  which  a fluctuation  of  either 


the  target  turn  rate  or  its  line  speed  is  noc  known  in  advance  (Figs.  4 and  6), 
the  step-by-step  control  action  is  taken  to  offset  the  estimation  error,  which 
prolongs  the  expected  intercept  intez^al.  Neveirtheless , reasonable  convergence 
occurs  after  a few  loops  of  estimation  in  these  particular  cases. 

Unknown  fluctuations  on  the  target  turn  rate  of  higher  an^litude,  e.g. 
-0.005+0. 005 (cos  0.2t+sin  0.2t)  wore  also  considered,  as  well  as  a time-varying 


Figure 

No. 

u^  (rad/sec) 

v^(ft/sec) 

2 

0.15 

1000 

3 

0.00 

1000 

4 

-0. 001+0. OOKcos  0.2t+sln  0.2 1) 

1000 

5 

0.00 

1000 

6 

0.00  980+20(cos  0.2t+sin  0.2t) 

sawtooth  target  line  speedy  in  order  to  test  the  effectiveness  of  the  astiM'* 
tlcn  sch«nM.  Sinulatlon  results  are  relatively  poor  In  these  cases  even  though 
further  estimation  steps  are  called  \;^on  in  trying' to  reduce  the  estimation 
error.  This  indicates  the  high  sensitivity  of  the  estimation  scheme  to  devia- 
tions fvaa  the  assumed  constant  values  for  the  unknown  parameters.  It  is  sug- 
gested that  a recursive  estimator,  which  utilizes  the  first  estimation  data  to 
initiate  a secondary  minimum  variance  estimation  on  either  the  target  tinti  rate 
or  line  speed,  may  be  worthy  of  consideration  for  futts^  analyses. 


VII.  Concluding  Remarks 

The  control  law  proposed  here  for  a missile  intercept  problem  is  admit- 
tedly ad  hoc.  Separating  the  estimation  and  control  problems,  then  combining 
their  solutions  in  a step-by-step  fashion,  does  not  in  any  way  constitute  an 


optimal  solution  for  the  over-all  problem.^ Moreover,  the  assumption  that  the 


pursuer  possesses  thriist  modulation  in  addition  to  the  usual  thrust  vectoring 
precludes  the  application  to  many  present  day  systems.  Nevertheless,  this 
assumption  does  allow  for  a single  closed  form  solution  to  the  minimum  control 
energy-terminal  constraint  problem  and,  it  is  believed,  has  enough  elements  of 
practicality  to  make  it  potentially  attractive  for  fut\n*e  systems  which  may 
possess  thrust  modulation  capabilities.  The  simulation  results  Indicate  good 
accuracy  for  an  intercept  when  the  estimation  errors  are  small.  At  the  same 
time,  the  simulations  for  various  amplitudes  of  the  fluctuations  in  the  target 
speed  and  turn  rate  reflect  the  possible  necessity  for  a higher  order,  or  more 
sophisticated,  estimation  scheme  to  alleviate  such  fluctuation  effects  on  termi- 
nal accuracy. 
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Figure  3.  Missile  Turn  Rate,  Speed  and  Relative  Trajectories 
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Figure  5.  Missile  T\jm  Rate,  Speed  and  Relative  Trajectories 
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Figure  6.  Missile  Tumi  Rate,  Speed  and  Relative  Trajectories 


|CORtTY  CLASSIFICATION  OF  THIS  PAGE  (When  l)titH^h.tttrrt*U)^ 

REPORT  DOCUMENTATION  PA6E  * 


n<  REPORT  NUMBER 


2.  GOVT  ACCESSION  NO. 


BSSM 


HEAD  INSTRUCTIONS 

DKI'ORE  COMPLF.TINO  FORM 

3 RECIPIENT’S  CATALOG  NUMBER 


AFOSR-TR-  78  - 0 0 2 2 , 


4.  TITLE  (mnd  Subllll*) 

CONTROL  LAW  FOR  A MISSILE  INTERCEPT  SYSTEM 


5.  TYPE  OF  REPORT  « PERIOD  COVERED 

Interim 


I 6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR!*; 

K.  C.  Wei  ' 
A . E . Pearson 


e.  CONTRACT  OR  GRANT  NUMBERr*} 

AFOSR  75-2793  ^ 


».  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Brown  University 
Division  of  Engineering 
Providence,  RI  02912 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Office  of  Scientific  Research/NM 
Bolling  AFB  DC  20332 


14.  MONITORING  AGENCY  NAME  « ADDRESSfl/  dlllmtutl  from  Conirelling  Olllea) 


10.  PROGRAM  element.  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 


61102F 

2304/Al 

12.  report  date 

Nov  77 

13.  NUMBER  OF  PACES 

23 

IS.  SECURITY  CLASS,  (ol  (hi*  nporl) 

UNCLASSIFIED 


IS*,  declassification/oowncrading 

SCHEDULE 


rieTDISTBIBUTrONSTAfEMENTTorihirRirirO 


Approved  for  public  release;  distribution  unlimited. 


• 7.  DISTRIBUTION  STATEMENT  (of  tht  abtlrmcl  mitrtd  In  Block  30,  II  dllltrtnl  from  Rtpotl) 


IS.  supplementary  NOTES 


fid.  KEY  WORDS  (Conilnu9  on  rovoroo  oido  i/  nocosoory  and  idmntify  by  btock  numbor) 


20.  abstract  fC^llnua  on  rovorao  aldo  1/  nocoaaafy  and  Idontliy  by  block  numbor) 

A feedback  control  law  is  derived  for  a two  dimensional  missile  intercept 
system  which  combines  in  a step-by-step  manner  the  closed  form  solutions  to  a 
least  squares  estimation  of  the  target  turn  rate,  line  speed  and  heading  with 
a minimum  control  energy  problem  subject  to  a terminal  intercept  constraint. 


OD  1 JAN*?!  1473  edition  OF  I NOV  SB  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  flWi*n  0*1*  Ifnlorod) 


ecUWITT  CcASSIFiC  A I .>J'.  Of  I nl&  PAOtfMTiwi  iJala  £nl«r*d; 


20.  ABSTlfA^X-(€^nh*tfed) 


The  major  assumptions  are  that  the  pursuer  possesses  thrust  modulation  capabi- 
lities, in  addition  to  thrust  vectoring,  and  that  continuous  measurements  of 

; 

the  range  coordinates  are  available  for  estimation  purposes.  Simulation  re- 
sults are  included  which  indicate  that  good  terminal  accurccy  is  achieved 
when  the  estimation  errors  are  small  and  illustrate  the  limitations  on  accu- 
racy when  fluctuations  occur  in  the  target  speed  and  turn  rate. 


unclassified 


SECUeiTy  CLASStFICATION  Or  PAGEfWh.n  Dmim 


